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Abstract 



We analyze the stationary problem for the Toda chain, and show that arising geometric data exactly 
correspond to the multi-support solutions of one-matrix model with a polynomial potential. For the 
first nontrivial examples the Hamiltonians and symplectic forms are calculated explicitly, and the 
consistency checks are performed. The corresponding quantum problem is formulated and some its 
If} ' properties and perspectives are discussed. 

o ' 



1 Introduction 

It was found long ago (see, for example [1, 2]) and emphasized recently [3, 4] in the context of string geometry 
and its certain consequences for the supersymmetric gauge theories, that the solution to one-matrix model 

Z = J* xJV W($)) (1) 



X 

5— ( 

■ with the polynomial potential 

W r (*) = ES*** fc (2) 

in planar approximation within the 1 /TV-expansion is geometrically described in terms of the hyperelliptic 
complex curve 

y 2 = W'(X) 2 — 4/ n _!(A) (3) 
endowed with generating differential (related to the eigenvalue density for the matrix <£>) 

dS = ydX (4) 

The parameters of potential (2) in this context are certain fixed constants or "Casimirs", while the dynamical 
variables are rather related with nontrivial periods of generating differential (4) on hyperelliptic curve (3). 

This geometric setup implies, in particular, that the most effective language for the planar solution of (1) is 
theory of integrable systems. Basically, already in [2] it was found that the partition function (1) T = ti 2 \ogZ 

1 Based on talks at "Classical and quantum integrable systems", Dubna, January 2005 and "Selected topics of modern math- 
ematical physics", St. Petersburg, June 2005, and a lecture for the minicourse: "Toda lattices: basics and perspectives", Fields 
Institute, Toronto, April 2005 
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at N — > oo with fixed HN is geometrically formulated as a (logarithm of) quasiclassical tau- function [5] . This is 
a particular case of the so called algebro-geometric " holomorphic" integrable systems, quite commonly arising 
recently in the context of string theory or Seiberg-Witten theory (see e.g. [11, 12] and references therein). 

On the other hand, it is also well-known that the matrix integral (1), as a function of parameters of the 
potential (2) and the size of the matrix N, is a tau- function of the (semi-infinite, forced) Toda chain [6]. The 
Toda chains are well-known families of integrable models arising in particular context of the algebro-geometric 
integrable systems: their afnne or periodic versions are literally the Seiberg-Witten integrable systems associated 
to the pure (TV = 2 SUSY) gauge theories [7]. One may be interested then, whether the integrable system, 
associated to (3), (4) has any relation to the Toda family. 

Below we are going to demonstrate, that the complex curve (3) with the differential (4) is directly associated 
with a particular integrable Toda system, namely with a stationary problem in the Toda chain 2 . We discuss the 
general structure of the classical stationary problem, present explicit computations for the simplest nontrivial 
examples where the "order n is n = 2 and n = 3, and discuss briefly the problems of its quantization. We hope 
that the observed equivalence of the matrix model and stationary Toda geometries will be first step towards 
understanding of proper analogs of the extra quasiclassical variables, and the considered below system can play 
the role of the corresponding "probe" model. 

2 Notations 

The Toda chain Lax operator is an infinite 3-diagonal matrix 

Lij = nSj^-i +p i 5 jl + r i+ iS jti+ i (5) 
(with i, j G Z), where it is convenient to introduce 

n = e |(9i-9i-i) (6) 

or 

R . = r 2 = gft-ft-i ( 7 ) 

to be useful in what follows. 

The Toda chain equations of motion follow from the Lax representation 

§i = [L,A 1 } (8) 

with Ai = TjlZ o L or, explicitly for the 7?.-matrix 

( A i)ij = \ (n+iS 3 ,t+i - nSj^) (9) 
providing (8) being equivalent to the well-known equations of motion 

for the chain of particles with the nearest neighbor exponential interaction. 

In a standard way, the nonlinear equations (8) can be formulated as consistency of auxiliary linear problem 

(Lip) i = r i+ iip i+1 + p t ip t + nipi-i = \ip l (11) 

2 Analogous stationary problem for the KdV hierarchy is well-known and widely discussed in the literature, see [8]. 
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The basis of two-dimensional space of solutions to the linear problem can be found in the form of expansion 

V± = A^e*** (!+<£ + £ + ...) (12) 
where for the coefficients of expansion one gets 

(13) 

C, - = Pi 

and 

(14) 

v'i -vr-i = Ri+i + & (£ -Ci) 

which relates them directly with the Toda chain phase space variables. 

3 Stationary problem 

An infinite Toda chain is a completely integrable system with (infinitely many) conserving charges, which 
generate the " higher flows" , commuting with (8) . The higher flows can be written as 

§- k =\L,A k ] (15) 

with the same TZ- matrix A k — \TZ o L k , as in (9). There are many ways to restrict the infinite Toda chain to 
a finite-dimensional subspace of its phase space; most well-known arc TV-periodic chain and open chain (the 
"Toda molecule") with N particles. 

Below we are going to concentrate on the stationary problem, determined by commutativity condition 

[L,A]=0 (16) 

with 

A = TTk=iCkA k (17) 
for some fixed number - "order" n, and some coefficients {cfc}. 

As any problem with commutativity of two operators, the stationary problem (16) acquires a complex curve 

dct (A(\) - Y) = (18) 

where matrix 

_ / o(A) 6(A) \ 
A{X) -{ c(A ) d (A) J 

with some polynomial matrix elements is representation of the operator (17) in the basis of two-dimensional 
space of the eigenfunctions (11) of the operator (5). The curve (18) is therefore hypcrelliptic 

y 2 = P 2n (A) = (a-df + Abc (20) 

with y = 2Y - TtA(X) = 2Y - a - d. 

The most natural symplectic form for the infinite Toda chain is just floe = X^ez*^ ^ ^Pi- ^ ' s a ^ so 
known that it "remains intact", when reduced to the 27V-dimensional phase space of the iV-periodic prob- 
lem or A^-particle open molecule, i.e. VLn — J^iLi^Qi A $Pi- The latter one can be also written as Qn = 
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Yes w=00 Tr 1 5L(w) A S^>)) in terms of the spectral parameter dependent Lax operator of the TV-periodic 
problem, see [9, 10, 11] for details. 

Similarly, the symplectic structure for the stationary problem can be determined by 

fi = res A =oo (dA 1* (* _1 <W(A) A**)) (21) 

where 6A(X) is now the variation (at constant spectral parameter A!) of the operator (17), (19) and 

* = ( 22 ) 

is the matrix of the Baker- Akhiezer functions, with being its matrix inverse. 

In what follows we will use for convenience a different from (12) normalization of the Baker-Akhiezer func- 
tions: 

/ \ 

(23) 

y o e — r- y 

Note, that the normalization of the Baker-Akhiezer function does not influence the symplectic form (21), since 
the difference, say for the change of normalization (23), is proportional to 

res (SYdX) A 5( qi - = (24) 

4 Explicit computations 

Let us now illustrate the general construction by explicit computations. 
4.1 n=l 

This is the trivial case of stationarity of the operator (9), which in the basis of eigcnfunctions (11) reads 

V V'i-i J \ n -±(\-pi-i) ) \ ipi-! ) 

giving rise to the genus zero curve (20) 

y 2 = A 2 - {pi + Pi _i)A + \{pi + Pt -i) 2 + 4R t (26) 

with the expansion 

ydX ~ rfA(A-|(p i +p i _i)+ 2 fi + ...) (27) 
The commutativity condition (16) [L, Ax] = leads to the further constraints 

Pi+i =Pi =Pi-i 
Ri+i = Ri 

and fixing the coefficients of (26) or (27) means that all variables turn into some constants. 
In the trivial n = 1 case, explicit calculation of (21) gives 

Ox = e^-^dqi a + \ (<*e+ a der + deti A rf d) - < + A d (29) 

Since 

i «+ a der + deti A d d) - < + a cii = I a d Pi + d Pi - x a (30) 

and due to (28) the symplectic form (29) vanishes. 



(28) 
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4.2 ri=2 

Consider now the simplest nontrivial case of cj, oc 5k n for n = 2, i.e. A = A 2 = \H o L 2 . Applying twice (5), 
one gets 

(i 2, 0), = r i+2 ri +1 ipi + 2 + r l+1 (p l+1 + p t ) fa +1 + (R i+1 + pf + Ri) fa + r, {pi + Pi-i) fa-i + nn-ifa-2 

(31) 

or 

(A 2 ip) i oc n +2 r i+ ifa +2 + r i+ i (p l+1 + p^ fa +1 - n (pi + Pi-i) fa-i - r i r l ^ 1 fa_ 2 (32) 

while stationarity of the gradient of the trace part TrL 2 = J2i {Pi + in (31), which disappears from (32), 
lead exactly to the constraints (28) of the previous trivial example. 

The commutativity constraints lead now to 

R i+1 (p l+ i +pi) = Riipi +Pi-i) 

(33) 

Ri+i +Pi = Ri-i +p1-i 

which reduce the dimension of the phase space to the four independent variables {Ri and pi on two neighboring 
sites) . 

Now, using (11), one can rewrite (32) in the basis of eigenfunctions of the Lax operator (5), e.g. 

/ * \ = fltX-p*-R i+ , + R i ) _ n{x + pi) w * \ 

\fa-i J V ri(A + Pi_i) -\{\ 2 -p 2 1 _ l -R l -i+R l ) ) \fa-i ) 

In particular, upon (33) one gets Tr^4 2 (A) = 0. 

(33) 

For the operator (34) the equation of the curve (20) reads 

y 2 = (A 2 -iQ) 2 -4i? 4 (A+ K )(A+^-i) 
Q = Ri+i + Ri-i - 2R t +pj+ p\_ x 
The expansion of the generating differential at A — > oo gives 

ydX x ~_ d\{\ 2 -\{Q + AR l )-2f{ Pl +p l - l )-^{2p l p l - 1 +Q + 2R t ) + ..) = 

d\(X 2 -C 2 -^-^ (piPi-i -Ri + C 2 ) - ^ - ^{C 2 + C 2 H) + ...) 



(35) 



A— *oo 



(36) 



Q+4R,=2C 2 



where fixing 2C 2 = Ri+i + R4-1 + 2R, t + p 2 + pf_ 1 together with the coefficient at ^ lead to the following 
Casimir constraints 

C'i = Ri(pi +Pi-i) 

(37) 

C* 2 = Ri+i +R l + p 2 = Ri-i+Ri+ p 2 -! 

which, in addition to (33) reduce the phase space to be two-dimensional. A convenient choice of co-ordinates is 
Ri (for some fixed site i) and the difference between the corresponding momenta p- = pi — Pi-\. 

The Hamiltonian H is, up to a numeric factor, a coefficient at the ^-term in (36) 

H = Ri {piPi-t -R t + C 2 ) = ? £ - \RiP- -Ri + C 2 R t (38) 

Pi-Pi-i=p- 4 "» 4 v ' 

With formulas (37) equation of the curve (35) acquires the form 

y 2 = (X 2 -C 2 ) 2 -4(C 1 X + H) (39) 
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and coincides exactly with the matrix model curve (3) for the potential W = 4p— C 2 X and function fi = C\X+H. 
Explicit calculation of the symplectic form for the n = 2 case, using (23) and (19), gives 

n 2 = \d (Ri-t + Ri + A d^ili - \d (Ri+i +R t + Pi) A d& + d( Pl - Pi _i) A dRi = 

= d(pi — Pi-i) A dRi ^ 

(37) 

Now, in order to fit between the hamiltonian (38), the curve (39) and the symplectic form (40) one should 
rewrite equation (39) in the " Seiberg-Witten" form (see e.g. [13] and [11]), i.e. as 

w + ^lA+ff = A 2 - C 2 

(41) 

w-^^- =y 
This means, in particular, that 2w = y + A 2 — C 2 , and 

\dy Ad\ = dw A d\ (42) 
The first of the equations (41) can now be written as 

H = w\ 2 - C\X - C 2 w - w 2 = w (A - <±f - % - C 2 w - w 2 (43) 
which exactly coincides with (38), provided by the identifications 

w <-> —Ri 



(44) 

„ ^ ^. — I— T» 

2 w 



A - & «- ±P 



4.3 n=3 

For the n = 3, with cx Sk3 the explicit expression for the operator A 3 = o L 3 looks like 

(A 3 ^) l oc r !+3 r I+2 r i+ i^ + 3 + r l+2 r l+ i(p l+2 + Pi+i + Pi)4>i+2+ 
{Pi+i + P1+1P1 +P 2 + r 2 +2 +R i+1 +Ri) Va+i - 

(45) 

-n (pj +PiPi-i +Pi_x + R4+1 + Ri + Ri-i) ipi-i- 
-ri-m- 2 {pi +Pi-i +Pi-2)4>i-2 - nn-in-^is 

while stationarity of the gradient of the trace part 

|ra 3 = §£ i p? + £ i p i (i2 i + i2 i+1 ) (46) 
gives rise to the constraints (33) of the previous example. The commutativity conditions (16) now give 
p\ + {Pi+i + 2^) Ri+i +p t R t = + (Pi_2 + 2pi_i) +Pi-iRi 

(47) 

^+1 + p%+m + Pi + R i+2 + Ri+i + Ri) = Ri (pI + PiPi-i + Pi-i + Ri+i + R l + R l -i) 

In the basis of eigenfunctions of the Lax operator matrix A3 (A) for (45) acquires the form 

/ ±(A 3 + 2r 2 A-p 3 -(p i+1 +2p i )r 2 +1 +p;_ 1 r 2 ) -r,(A 2 + Pl X + p 2 + r 2 +1 + r 2 ) \ 

\ n(X 2 +p l _ 1 X+p 2 _ 1 +r 2 +r 2 _ x ) -±(A 3 + 2r 2 A - p\_ x - ( Pl _ 2 + 2p l _ 1 )r 2 _ 1 + p t r 2 ) J 

(48) 
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For the operator (48) the curve (20) reads 

y 2 = A 6 - (p? +p?_! + (pi_2 + 2p,_ 1 )i? 4 _i + 3{ Pt + Pl -i)R t + {p i+1 + 2 Pi )R i+1 ) A 3 - 
-4ifc + i? 4 + + p 2 + p^i-i + A 2 - 2R t ((p i+1 + 2 Pl + 2p l _ 1 )R l+1 + 

+(Pi+Pi-i)Ri + {pi-2 + 2pi-! + 2 Pi )R l - 1 + p 3 + 2p\p^ 1 + 2p i p?_ 1 +pti) A (49) 
+3 ((Pi+i + 2pi)R i+1 - {pi+pi-^Ri + (pi-2 + 2p i _ 1 )R l - 1 + pf +pf_i) 2 - 
-4i2< (p 2 + + «i) (p 2 _! + iZj + Ri-i) 
The expansion of the generating differential is 

yd\ ~ dX (A 3 - | ((p i+ i + 2 Pl )R l+1 + 3(p, + p t ~ 1 )R l + ( Pl - 2 + 2 Pl - 1 )R t _ 1 + pf + p 3 _ x ) - 

+ i?« + J^_! + p 2 + pip^-L +pj_ 1 ) - %t ((Pi+i + 2 Pl + 2p i - 1 )R i+1 + 
+(Pi+Pi-i)Ri + {Pi-2 + 2pi_i + 2pi)Ri-i +P 3 i+ 2p\ Vl - X + 2p i p?_ 1 +pf_ 1 ) - 
-ft- ((Pi+iPi + Pi+iPi-i + 2 P f + 2piPi_ x + 2p 2 _ 1 )R l+1 + (3p 2 + 3p 2 _! + 2 PiPl - 1 )R i + 
+{Pi-iPi- 2 + PiPi-2 + 2 PiPt - 1 + 2p 2 _ 1 + 2p 2 )R i _ 1 + 2{R l+1 Ri + R4+1R4-1 + B 2 + RiRi-i)+ 
+pf + plpi-! + 2p 2 p 2 _ 1 + p ;i p\_ x +pf_ 1 ) + ...) 

Fixing the coefficients at singular terms of the expansion of ydX to be the independent upon dynamical variables 
numbers (or "Casimirs"), using (48), one can write 

Ci = Ri (p 2 + PiPi-i + p 2 _x + R i+ i + R t + Ri-i) 

(51) 

C 3 = Pi + (Pi+i + 2 Pi ) R i+ i + (pi-i + 2 Pi ) R, = p 3 _! + (p»_ 2 + 2p;_i) + (2p,_! + p^ R, 

Together with (47), relations (51) impose four constraints to the eight-dimensional space of variables (Rj and pj 
with j = i — 1, i, i + 1 together with Ri-2 and Pi+2), so that the phase space in this example is four-dimensional. 

Then the expansion (50) upon (51) turns into 



(50) 



ydX ~ dX (A 3 - C 3 - ^f 1 {Ri+i + R t + R4-1 +pj+ P*P*-i + P-_i) - 

A—* 00 

(pi-iRi+i - {pi+Pi-^Ri+PiRi-! + P 2 p t -i +PiP 2 _i + C 3 ) + ...)= (52) 
= rfA(A 3 -C 3 -^-^-^ + ...) 

so that the first of two independent hamiltonians (proportional to the coefficients at |£ and j^) can be, using 
(51), presented as 

Hi = C 3 R, + %-ipi+Pi-t) - \Ri (p!+pt 1 + 3R l (p t +p l -i) + (Pi - Pi-i)(Ri+i - Ri-i)) (53) 

and 

H ^ = jk + ^t - Pift-i) + C 3 RiijPi + Pi-i) + ±i? 3 - 

-\Ri{R i+ i - i?,-i) 2 - \Ri{Ri+i - R t -i)( P 2 - P ti)- ( 54 ) 
-\Ri {pj + 3p?p?_i + pU) - Rt (p? + IPiPi-i + Pti) 
The equation of the curve (49) now becomes 

y 2 = (X 3 -C 3 ) 2 -4(C 1 X 2 + H 1 X + H 2 ) (55) 
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again coinciding with (3) with the potential W 



C 3 X and function / 2 = CiA 2 + #iA + H 2 . 



Explicit calculation of the symplectic form gives 

^3 = ^C-l A d {Pi-1 + (Pi-2 + 2Pi-l)Ri-l + {Pi + 2p i _ 1 )R i ) + 

+ \d (pf + ( Pl+1 + 2 Pi )R i+1 + (pi_! + 2pi)Ri) A dg 
+d(R i+ i - B4-1) A dR, + 2Ridpi A + 
+rf (p? - A + (jpi-idpi - p l dp i - 1 ) A 

Using (51), and introducing i?_ = iij+i — i£,-i this slightly simplifies to 

CI3 = A dRi + 2Ridpi A c?p s -i + 
+rf (p? - A dR t + (pi-xdpi - Pidpi-{) A dR t 

(certainly, d(Cls) = 0), or 



f^3 = (dR- dRi dpi dpi-i) ■ ^3 



/ dR- \ 
dRi 
dpi 
V dp^i J 



where the wedge product is implied and ^3 is the matrix 

/ \ 

~\ -(^ + £ f i ) ft-i + f 

p 4 + £ r i ^ 

V -( Pi _i + f) -ifc 

A nontrivial check is that the Hamiltonians (53) and (54) indeed commute 

{Hi,H 2 }n 3 = 



(56) 



(57) 



(58) 



(59) 



(60) 



with respect to the Poisson bracket, corresponding to symplectic structure (57) or defined by the inverse to (59) 
matrix 

/ 



V 

giving rise to the "elementary brackets" 






-2 


2pi-i+Pi 

Ri 


2pi+pi-x 

R, 


\ 


2 













2pi-i+p; 

Ri 








1 

Ri 




2pi+p s -i 

Ri 





J_ 

Ri 





/ 



(61) 



{Ri, R-}n 3 — 1 
{Pi,Pi-i}n 3 = -2R- 
{R-,Pi}n 3 = {Pi-i + f ) 
{R-, Pi -i}n 3 = w i {pi + Ei 2 zl ) 



(62) 
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The symplectic structure (57) can be brought to canonical form, introducing new variables as pi = p cosh# 
and pi_ i =p sinhfl, then 

fl 3 = dR- A dR, + Rid(p 2 ) A d9 + p 2 dR, L A dO + d(p 2 ) A dR, = 

(63) 

= d(i?_ + p 2 ) A dRi + d{R iP 2 ) A dO 
and one concludes that the Darboux co-ordinates are 6, R = Ri together with 

P = R*P 2 , N 

64 

A = R_+p 2 

In these variables the Poisson brackets are 

{A,i?K = l, {p,6}n 3 = l (65) 
while all other vanish. The Hamiltonians (53), (54) in these variables have the form 

Hl = C 3 R + Q±p^R-V 2 e e - lApLpRWe-* + § (e- e - e 36 ) p^R- 1 ' 2 - |p 1 /2 i? 3/2^ 
H * = m + f R + f i ( e ~ W - e29 ) + Csp^R^e 6 + ii? 3 - \A 2 R l P Re 29 + (66) 

+i pR e-*> + ±£(e«> + e-«>) 

i.e. are generally functions of the fractional powers of dynamical variables. The dependence on fractional 
powers, however, disappears for H 2 at vanishing Casimir C3 = and for the square H 2 at C\ = C3 = 0. 

5 Discussion 

Let us, finally, discuss some related issues and open problems. 

Geometry on stationary and periodic problems 

When defining the stationary problem, we have started with an infinite Toda chain. Remaining all ingredients 
almost intact, one could take as initial point a TV-periodic problem in Toda instead (qi+N = Qi, Pi+N = Pi, 
tpi+N = wipi etc) with sufficiently large N. 

Such procedure, however, may be already treated purely in terms of algebro-geometric integrable systems. 
One may start with the curve of TV-periodic chain 

W + ± = P N (\) (67) 

endowed with generating differential dS — A^Jr an d consider its reduction down to 2n-dimensional subspace 
in the moduli space for n < N. A natural way to impose such reduction is to require existence on the curve 
(67) of a single-valued meromorphic function with only two poles (at A = 00) of order n. Existence of such 
function on the curve (67) is exactly equivalent to stationarity of the order n flow in periodic Toda, associated 
with the linear combination of meromorphic differentials dfiA = Sfc=i c kd£lk- Since there is no natural function 
with desired properties on the curve (67) in general position, this constraint effectively reduces the (smooth) 
genus of the Ricmann surface from N — 1 to n — 1 and one ends up with the curve of a stationary Toda chain 
(3). Such reduction was already discussed in the literature comparing the Scibcrg-Wittcn and Dijkgraaf-Vafa 
geometries, see [3] and, for example, [14, 15, 16]. Note also, that the corresponding generating differentials or 
symplectic structures remain unrelated by this procedure, and this exactly corresponds to the difference between 
the symplectic structures (21), (40) and (57) of stationary problem and that of the infinite chain we discussed 
above. 
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Remarks on quantum case 

For the curve (39) and the symplectic form (42) one can immediately write the "naive" Schrodingcr equation 

h2 & T = (W{\) 2 - hW"{\)) T (68) 

with the solution 

T = exp(-IW(A)) (69) 

since 

h2 & W'(\) 2 + hW"(X) = {h£ W'(\)) {hJL + W'(\)) (70) 

Comparison of (68) with (39) for the potential W'{\) = A 2 — C2 gives H = and C\ = l|, and comparison 
with (55) leads to Hi = H\ — and Ci = Formulas (41) and (42) rather suggest that the quantization 
can be better performed in the (y, u>) -variables. For the periodic Toda chain (67) this way leads to the Baxter 
second-order difference equation, but for the stationary problem equation (41), due to (42) would give rise to 
an integral operator. The main obstacle on this way is that it is vague what are the "proper" variables for 
quantization of stationary problem, since using the naive Darbough co-ordinates, already in n — 3 case, leads 
to the non-analytic dependencies in the Hamiltonians. 

One can nevertheless try to impose the quantum relations 

H l Z = H l Z, i=l,...,n-l (71) 

onto the matrix model ^-function, which instead of naive partition function Z defined by integral representation 
(1) contains explicitly information about (quantum version of) geometry (3), (4), i.e. one literally has to solve 
(71) for the function 

Z{c,H) = exp(^F {c,H) + ...) (72) 

This formulation of quantum problem is consistent with the fact that the form of Hamiltonians (66) simplifies 
at vanishing (clasically) Casimirs, when they are of order of h or being some quantum corrections. The classical 
formulas (66) for Hi (or some polynomial functions of them) are then replaced by differential and difference 
operators Tti via A — ► d/dR and 9 — ► d/dp, so that e ±e turn into the shift operators in p- variable. Constructed 
in this way Z-function (72) can be treated as a matrix-model analogue of Nekrasov's generalization of the 
Seiberg-Witten prepotential [17]. 
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